:00 


- — < — 

.00 

0                 1Г) 

?: 1Г) 

х=^0 

0=^00 

:.    О 

-    1 

^ о 

т 

г СО 

^= [>- 

1 

СО 

ЫиЫтоуа,   Е  А 

О  гас11оак-Ыупот  гагоёгеуе 
ктрпукЬ  аБ1его1йоу 


ОБ 

153 


два  А 


А(Зт^%)ЛШЧЕСКиИ  ЖУРНАЛ 


^^  Е.  А.  ЛЮБИМОВА  и  А.  Г.  СТАРКОВА- 

о  РАДИОАКТИВНОМ  РАЗОГРЕВЕ  КРУПНЫХ  АСТЕРОИДОВ 
И  СТРУКТУРЕ   МЕТЕОРИТОВ 


Приводятся  результаты  расчетов  изменения  температуры  в  центре  тел  разме- 
рами крупных  астероидов  (о?=  100  —  400  км),  обусловленной  теплом  радпоактпв- 
ного  происхождения,  в  предположении  равного  содержания  радиоактивных  элемен- 
тов в  земном  и  метеоритном  веществе. 

Делается    вывод    о   том,    что    в   далеком  прошлом   (3 — 5.10'   лет   тому   назад) 

V.   таких   телах    могло   накапливаться   достаточное   количество  радиогенного  тепла, 

пособного  разогреть  их  вещество  до  высокой   температуры,    и  что,  таким  образом, 

южно  объяснить  следы  разогрева   и   расплавления   в    метеорит ^>      -.г-п?  счптпть  их 

сколками  подобных  тел. 

Космогоническая  теория  академика  О.  Ю.  Шмидта  [1]  утверждает  на- 
личие единого  процесса  происхождения  и  развития  для  всех  тел  солнеч- 
ной системы:  планет,  астероидов,  метеоритов  и.  комет.  Исследование  про- 
цесса образования  планет  из  газо-пылевого  облака  [2]  приводит  к  заключс- 
!1пю,  что  сперва  образуется  много  небольших  тел  с  массами  порядка  масс 
I стероидов.  В  связи  с  этим  высказывалась  гипотеза  [3,  4],  что  метеориты 
являются  осколками  подобных  тел,  не  входивших  в  состав  какой-либо 
крупно11  планеты.  Однако  структура  метеоритов  (особенно,  железных) 
носит  на  себе  следы  разогрева  и  даже  расплавления.  Поэтому  возникает 
вопрос,  могли  ли  тела  размером  с  астероид  разогреваться  до  пыгокоГг 
температуры. 

Как    показано    в    статье    [5],  незначительная    примесь    радиии1х1тяп.1А 
эле.ментов,  обнаруженная  в  зе.мных   породах  и  метеоритном  веществе,  обе- 
течивает  сильный  разогрев  Зе.мли  за    3 — б   хмлрд.    лет  ее  существования. 
Поскольку  радиоактивные  элементы  непрерывно  распадаются,  то  несколько 
миллиардов    лет    тому    назад,    в    эпоху    образования    планет,    выделение 
[радиогенного  тепла  было  в  несколько  раз  больше,    чем   теперь.  Величина 
первоначального  содержания    радиоактивных   элементов,    вычисленная    на 
основании   данных   об    их   современном   содержании,    зависит  от  возраста 
тела.    Че.м  больше    возраст,   тем    больше   их   первоначальное    содержание. 
"    По  расчетам  Рамзея  [6]  и  Буллена  [7)  планеты  земной  группы  имеют, 
основном,  одинаковый  состав,  совпадающий  или  почти  совпадающий  со 
дним  составом  метеоритов.    Сходство   изотопного   состава   ряда  веществ 
1МЛИ  и  метеоритов  указывает  на    то,    что    они  образовались  из  одного  и 
(ГО  же  первичного  вещества.  Это  дает  основание    считать,  что  и  средшч' 
(Держание  радиоактивных  элементов  в  Земле  не    слишком  отличается  от 
1'теоритного.  Такое  же  содержание  должно   быть  в  других  телах  земно11 
1>уппы,  к  которым  следует  отнести  и   астероиды.    По,  хотя  это  содержа- 
ние было  достаточным  для  разогрева  Земли,  без  специального   рассмот|)е- 
пия  нельзя  установить,  достаточно   ли   оно   для  разогрева  малых  тел,  из 
которых  тепло  легче  ускользает  в  пространство. 

С  этой  целью  нами  были  пр^зделаны    расчеты    температурного   режима 
1ЛЯ  шаров  радиусом  от  200  до  50  км  (размеры  крупнейших  астероидов), 
-"•условленного    радиоактивным    разогревом.    Предполагалось,    что    радио- 
;ич'тивные  элементы  равномерно  рассеяны  в  этих  телах,  а  плотность,  тепло- 
проводность и  теплоемкость  всюду  одинаковы  и  постоянны. 
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Для  тел  малой  л'ассы  эти  предположения  справедливы,  так  как  (если 
не  считать  мелкие  неид породности)  вещество  этих  тел  вряд  ли  могло  под- 
вергнуться значительному  рассл1>ению,  а  сравнительно  небольшое  давле- 
ние не  могло  СИЛ1  но  изменять  свойства  вещества  внутри  этих  тел. 

Температура  поверхности  полагалась  постоянной.  Как  показывает  про- 
стой расчет,  она  не  могла  заметно  изменяться  вследствие  разогрева  от 
внутренних  источников  тепла. 

Действительно,  если  бы  все  тепло,  выделяемое  внутри,  излучалось 
в  пространство,  то  температура  поверхности  определялась  бы  выражением 

Т^у   -7^,  где  а — константа  Боль:(^ана,  равная  5.72-10~^эрг  см"2  сек.~^ 

град."^,  К  —  радиус  шара,  ^  —  количество  тепла,  выделяемое  в  1  см^  в  1  сек. 
Полагая,  что  в  настоящее  время  ^  —  ^^  ~  50-10"^^  кал/см^  сек.,  а6-10^лет 
назад  д  =  д2  =  ^ЬО'И)^^^  кал/см^  сек.  (см.  фиг.  1  в  статье  [5]), 
получим  при  Я  =  200  км,  2'(^1)  :==;  20°К,  /"(дз);^  40°К.  На  самом  деле  Т 
будет  еще  меньше,  так  как  в  пространство  излучается  не  все  выделяемое 
внутри  тепло.  Следовательно,  поверхностная  температура  астероидов,  как 
и  Земли,  определяется,  в  основном,  действием  Солнца,  не  проникающим 
глубоко  внутр]>  и  распределяющимся,  в  среднем,  равномерно  на  поверх- 
ности тела  вследствие  его  вращения  вокруг  оси. 

Начальная  температура  полагалась  равной  нулю.  С  точки  зрения  космо- 
гонической схемы  образования  тел  солнечной  системы  путем  объединения 
холодных  частиц,  это  предположение  вполне  справедливо  для  тел  малой 
лтссы,  так  как  интенсивность  процессов,  обуславливавших  первоначальную 
температуру  (удары  частиц  о  поверхность  и  сжатие  внутренних  частей 
за  счет  роста  внешних)  для  таких  тел  незначительна.  '1аким  образом, 
постановка  задачи  остается  такой  же,  как  в  статье  [5],  и  решение  ее 
представляется  в  виде: 

О  о 

2  [ф(^11±М1^^)^ф({^ШМ^)]а.,  (1) 

;Йо^    V  У/./с(/-т)  у        V /4Л(г-т)  у/    '  ^  ^ 

где  Т  —  температура,  г  —  расстояние  от  центра  тела,  ^  —  время  с  момента 
начала  разогрева,  к  —  температуропроводность,  с  —  теплоемкость,  р  —  плот- 

—  г—  2  *  1 

пссть,  1{  —  радиус  тела,  Ф  {х)  —  интеграл  ошибок     Ф{х)—  1  —  -уг-  Г  е-^'с^а 

о 

т 

И{1)  —  ^Н1е~^  — количество    тепла,    выделяемого   в   единицу    времени 

г-1 

и  единице  обьема  всеми  радиоактивными  элементами,  для  каждого  из 
которых  ).^  —  константа  распада,  //{  —  первоначальное   содержание. 

Формула  (1),  примененная  в  [5]  для  случая  /?  =  Язсмлн,  вполне  при- 
ложима  для  расчета  температуры  ьнутри  сфер  меньших  радиусов,  так 
1хак  ряд,  стоящий  под  интегралом,  достаточно  быстро  сходится  и  для  этих 
случаев.  Дейс1Вительно,  функция  Ф  (дг)  является  очень  быстро  убывающей 
(  ростом  аргумента:  при  х—'б,  Ф  (3)  имеет  порядок  малости  10"^,  при 
•/•  =  4,  Ф(4)— 10"^.  Так  как  множители,  стоящие  перед  Ф{х),  имеют 
порядок  —10^,  то  для  вычисления  Т{г,1)  достаточно  ограничиться  чле- 
нами, содержащими  Ф(а;),  при  а;<3  и  отбросить  члены  с  л;>-3.  Аргу- 
^^ентом  х  является  функция 

х=:г„=<А"+Дг±-'.  (2), 
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11'  I  м!,)  ),  1  I  \лен1шснием  /?  для  определения  Т  вблизи  центра 
(/•--О)  нужно  брать  Бсе  бильшее  число  членов.  Денствител! но,  положив 
т  =  0  и  фиксируя  наибольшее  значение  ^  =  710^  лет  (рассматриваем  / 
п  достаточно  большом  интервале  0</<:7  108  лет),  найдем  наибольшее 
еловое  значение  знаменателя  в  выражении  (2)  при  Л-  =  0.01 


УАк  {{ —  т)  =  /4 •  0.001  •  7  •  ;11  Г)  •  1  < )>«  =  V 0«  к-м. 

/9;.        I        \\      К 

неравенства  Хп=  ^~   ^.^  — ;>  3    следует,    что    если    /?.>2700    км,   то 

можно  ограничиться  одним  первым  членом  формулы  (1);  если  /?-^Г00км, 
го  можно  отбросить  все  члены  ряда,  за  исключением  первых  двух;  если 
/^--'ЮО  км,  то  следует  принимать  во  внимание  члены  до  п=  13.  Указан- 
ные оценки  сделаны  в  предположении,  что  коэффициент  температурс'про- 
водностп  /и  =  0.01.  При  меньших  значениях  к  можно  ограничиться  еще 
меньшим  числом  членов.  Так,  например,  при  к  =  0.002  и  Л  =  100  км 
необходимо  принимать  во  внимание  члены  лишь  до  /г  =  5. 

Рассмотрение  предельного  перехода  в  формуле  (1)  при  г->0  приводит 
к  более  простой  формуле  для  вычисления  температуры  в  центре  тел: 

,  со       -  /<2»+1'Д  \  2 


=0 


пли 


=0 


I 


Однако  число  вычислительных  операций  при  использовании  формулы  (3) 
практически  не  уменьшается.  Поэтому  безразлично,  по  какой  из  формул  (1) 
пли  (3)  вести  расчеты. 

Прнводи.мые  ниже  результаты  получены  путем  численного  интегриро- 
наиия  выражения  (1)  при  малом  значении  г. 

Содержание  радиоактивных  элементов  бралось  при  расчете  равным 
одной  из  оценок  среднего  содержания  их  в  метеоритном  веществе  (кон- 
центрация Б  в  [5]): 

Уран     Ы0-'  г/г. 
Торпй  3.00- 10-'  г/г. 
Калпй  0.15- 10"^  г/г. 

(При  концентрации  А  приводимые  ниже  температуры  нужно  уменьшить 
примерно  в  2  раза). 

Константы  распада  и  тепловыделения  /.{  и  //^  взяты  здесь  такими  же, 
как  в  [5].  Телшературопроводность  к  для  малых  космических  тел  должна 
быть  меньше,  чем  для  Земли  (в  [5]  к  полагалось  равным  0.01  С.С8),  так 
Вк  вещество  этих  тел  почти  не  сжато  и  является,  вероятно,  бо.тее  рых- 
Шш,  чем  поверхностные  породы  Земли,  для  которых  к  =  0.00(5  :|1.>(  I.  /.• 
принято  рапным  0.002. 

Согласно  последным  определениям  Рик  [8]  по   данным  аргоиоиою  ли- 
ла,  вещество    метеоритов    существует    в    твердом  состоянии,  по  крайней 
,.ре    4 — 4.5-10^  лет.    Полный    возраст    метеоритов   с    учетом    начального 
гапа,  когда  вещество  могло  пройти  стадию  значительных  температу1>иых 
изменений,  должен  быть  больше.  В    расчетах    возраст   /  принимался  рав- 
]1ым  4  и  6-10^  лет. 

Па  рис.  1  изображено  изменение  температуры  со  временем  в  центре 
тела  радиусом  200  км,  полученное  в  предположении  разного  возраста  /. 
На  рис.  2    /  =  6-10»  лет,  а  Л  =  50,  100  и  200  км. 

Рис.  1   и  2  показывают,  что  на  ранней  стадии  существования  рассмат- 
ваемые  тела  претерпели  сп.'г.пьпЧ    разогрев.    Температура  в  нехитро  .'•тп\- 
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тел  достигала  максимума  в  первые  полтора  миллиарда  лет.  К  настоящему 
времени  эти  тела  являются  остывшими. 


5000- 


4  I  в  10" лет 


Рис.  1.     Изменение      температуры      со 

временем  в  центре  тела  радиусом  В.  = 

=  200  км;    к  =  0.002,  концентрация  Б; 

^_^_6.10^»;  2  — /  =  4-10» 


О        1 
(ВЮ^лт  тщ 

Рис.  2.  Изменение  температуры  со 
временем  в  центре  малых  тел  разных 
Л;А:  =  0.002;  /^б-Ю»;  концентрация[Б; 
7  — Д  =  200  км;  2  — Л  =  100  км; 
5  —  Л  =  50  км 


В  табл.  1  приведены  максимальные  значения  центральной  температуры 
Гщах  и  моменты  времени  1^,  в  которые  она  достигалась  при  разных  1яН. 


I 

я  п  км 

6-10» 

лет 

4.10* 

лет 

Т'  тах 

г,.1Г' 

Т"  тах 

<,.10-» 

200 

100 

50 

4800 

1700 

400 

1.1 
0.4 
0.3 

1600 
500 
130 

1.2 
0.5 
0.3 

Проделанные  расчеты  показывают,  что  в  далеком  прошлом  в  наиболее 
крупных  телах  астероидного  типа  (радиусом  от  50  до  200  км)  могло  на- 
капливаться достаточное  количество  радиогенного  тепла,  способного  разо- 
греть их  до  высокой  температуры.  Таким  образом,  гипотеза  о  том,  что 
метеориты — это  осколки  тел  астероидных  размеров,  никогда  не  входивших 
в  состав  крупной  планеты,  не  противоречит  признакам  разогрева  и  расплав- 
ления, обнаруженным  в  структуре  метеоритов. 

Авторы  выражают  благодарность  Б.  И.  Давыдову  и  Б.  Ю.  Левину 
за  ценные  замечания. 
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